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Abstract 

We report a high field investigation (up to 45 T) of the metamagnetic transition in Celrlns 
with resistivity and de-Haas-van-Alphen (dHvA) effect measurements in the temperature range 
0.03-1 K. As the magnetic field is increased the resistivity increases, reaches a maximum at the 
metamagnetic critical field, and falls precipitously for fields just above the transition, while the 
amplitude of all measurable dHvA frequencies are significantly attenuated near the metamagnetic 
critical field. However, the dHvA frequencies and cyclotron masses are not substantially altered by 
the transition. In the low field state, the resistivity is observed to increase toward low temperatures 
in a singular fashion, a behavior that is rapidly suppressed above the transition. Instead, in the 
high field state, the resistivity monotonically increases with temperature with a dependence that 
is more singular than the iconic Fermi-liquid, temperature-squared, behavior. Both the damping 
of the dHvA amplitudes and the increased resistivity near the metamagnetic critical field indicate 
an increased scattering rate for charge carriers consistent with critical fluctuation scattering in 
proximity to a phase transition. The dHvA amplitudes do not uniformly recover above the critical 
field, with some hole-like orbits being entirely suppressed at high fields. These changes, taken as a 
whole, suggest that the metamagnetic transition in Celrlns is associated with the polarization and 
localization of the heaviest of quasiparticles on the hole-like Fermi surface. 
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Itinerant Electron Metamagnetism (IEM) refers to a field induced transition from a para- 
magnetic (PM) to a field-polarized paramagnetic state, first predicted for exchange-enhanced 
paramagnets such as Pd^. The effect is intriguing since large, almost discontinuous changes 
in the magnetization can be observed through such a transition. In fact such field-induced 
magnetization jumps have usually been interpreted as a dramatic change in the density of 
states (DOS) at the Fermi energy as a Zeeman split DOS peak crosses the Fermi level. 
The most spectacular examples of IEM have been observed in the MgCu2 type cubic Laves 
phases such as AC02 (^4=Sc, Y, Lu)-. In heavy fermion systems, the screening of local 
moments by conduction electrons leads to a Kondo resonance peak in the DOS close to the 
Fermi level, so that large magnetic fields can lead to IEM behavior. Among the first heavy 
fermion compounds reported to have a metamagnetic transition are UCoAl^, and UPt^, 
and, by far the most thoroughly investigated, CePa^Siz 5 . Theoretical approaches that take 
into account the spin fluctuations can successfully capture the basic features of IEM in Laves 
phase compounds^. In order to capture some of the more subtle aspects of IEM in heavy 
fermion systems, Hubbard or Anderson model physics must be included. These models can 
reproduce features such as the Fermi surface volume change inferred from the dHvA- and 
due to /-electron localization 8 , as well as the sign change of the exchange coupling 9 (from 
antiferro- to ferromagnetic) observed in neutron scattering^. 

Despite a great amount of both theoretical and experimental work, some of the essential 
questions about IEM remain unsolved. For instance, is there a common picture emerging 
between the d- and /-electron systems that exhibit IEM? Would it involve solely a change 
in Fermi surface geometry (along with a magnetoelastic instability) or include a change in 
the volume enclosed by the Fermi surface? This is a particularly relevant issue for the heavy 
fermion systems, since Fermi surface changes have been observed under applied fielcUi^ and 
pressure^ and it is not clear if this indicates /-electron localization due to field polariza- 
tion of the conducting electrons. Finally the discovery of IEM in the Ruddlesden-Popper 
compounds Sr 3 Ru207 and Sr 4 Ru30io has opened new horizons^— . Systematic deviations 
from the Fermi Liquid (FL) theory are reported in the resistivity (p)^, thermal expansion^, 
and NMR— near the metamagnetic transition in Sr 3 Ru207, raising the possibility of a field- 
induced quantum critical point (QCP) in IEM systems. The T = suppression of a first 
order metamagnetic discontinuity terminating at a critical end-point has been shown the- 
oretically to result in a QCP—, a situation potentially realized in Sr 3 Ru207. The role of 
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disorder— >2i as well as the formation of magnetic domains^ in the proximity of a metamag- 
netic transition are still under investigation. 

Here we report a high field investigation of the recently discovered IEM transition in 
Celrln^ - — , to search for evidence of a metamagnetic QCP. Celrln 5 is a heavy fermion 
superconducting member of the layered tetragonal CeMIn 5 (M=Rh,Ir,Co) family-. Since 
their discovery, the 1-1-5 compounds have been the focus of much attention because they 
exhibit striking deviations from FL theory. There has been a great deal of work trying 
to understand the origin of such anomalous behavior in strongly correlated systems^ and 
descriptions are often based on the proximity to a QCP in their phase diagrams. Near 
such a singular point where a magnetic phase transition is driven to T = by an external 
tuning parameter (pressure, magnetic field or doping) the presence of quantum fluctuations 
is expected to have a strong effect on all physical properties, leading to the breakdown of the 
FL theory. In particular, the 1-1-5 compounds have rich and complex phase diagrams, with 
notably a field tuned QCP near the superconducting upper critical field in CeCoIn^ - — and 
a pressure tuned QCP in CeRhln^^i Besides their phase diagram, their heavy fermion 
state has also attracted attention in its own right. A careful analysis of the transport and 
thermodynamic properties of La-diluted CeCoIn 5 and Celrln 5 suggests that the coherent 
heavy fermion ground state coexists with a finite fraction of single-ion Kondo centers down to 
the lowest T— . This conclusion was further supported by optical conductivity measurements 
on CeMIns samples that revealed a broad hybridization between the conduction electrons 
and the more localized /-electrons which included contributions from 4 electronic bands^. In 
addition, the hybridization gaps that result appear to have extensive momentum dependence 
with regions where the smallest gap, associated with a particular hole band ((3 orbits in de 
Haas- van Alphen experiments), is seen to go to zero. The conclusion was that Kondo 
screening could be incomplete, particularly in regions where the hybridization gap goes 
to zero. Band structure calculations using DMFT, motivated by the optical conductivity 
measurements, have found two distinct hybridization gaps in Celrlns, attributed to two 
distinct In sub-bands hybridizing with the Ce lattice^. 

In addition to the many common trends with its Co and Rh counterparts, a distinguish- 
ing feature of Celrln 5 is the presence of a metamagnetic-like transition near 30 T for the 
[001] field orientation (perpendicular to the Celn3 planes). This transition is characterized 
by a non-linear increase in magnetization^^ and a A-like anomaly in the specific heat^. 
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Because the family of 1-1-5 compounds has many unusual features, including what appear 
to be unusual transitions perhaps associated with interesting QCPs, we are compelled to 
explore the possible existence of a metamagnetic QCP in Celrln 5 . Our earlier specific heat 
and p measurements up to 17 T have shown a field-induced non-Fermi Liquid behavior 
(NFL). However the limited field range of these experiments did not allow us to explore the 
immediate vicinity of the metamagnetic transition^. We have subsequently extended our 
investigations to higher fields (up to 45 T) and explored the angle dependence of the high 
field transition in Celrlns combining resistivity and torque magnetometry measurements^. 
Here, we give a complete account of our high field investigations with an in-depth analysis 
in search for Fermi surface changes at the transition and for clues about the cause of the 
NFL behavior. 

Our measurements reveal several changes that occur at or near the metamagnetic critical 
field which have important consequences for interpreting the mechanism for IEM in Celrln 5 . 
The most obvious is that an increased scattering rate for charge carriers results in an in- 
creased resistivity and a decreased dHvA amplitude in the region of the transition. This is 
accompanied by a resistivity anomaly at very low temperature that switches from increasing 
resistivity with decreasing temperature below the transition, to decreasing resistivity with 
decreasing temperature above the transition, none of which can be interpreted as a Fermi 
liquid-like behavior. Finally, we observe the complete suppression of two of the heaviest 
dHvA orbits in the high field state, while other orbits recover and eventually exceed their 
lower field amplitude. These changes are not accompanied by significant variations in either 
the dHvA frequency, indicating no measurable change to the extremal areas of the Fermi 
surface, or the cyclotron masses associated with orbits that we observe. Our data indicate 
that light parts of the Fermi surface are left unchanged by the IEM while the heavier parts 
are suppressed beyond the metamagnetic critical field. Although our data do not provide 
direct evidence for a metamagnetic QCP, we observe NFL transport properties in proximity 
to the critical field as well as an instability of the heaviest Fermi surface sheet, both of which 
are consistent with quantum critical behavior. 

The paper is organized as follows: we first present the high field phase diagram of Celrln 5 , 
which has been previously established via magnetization^ and specific heat measurements^. 
This is followed in section III by a description of the resistivity anomalies associated with 
the metamagnetic transition. Section IV is devoted to the evolution of the Fermi surface 
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across the metamagnetic transition, with an analysis of the dHvA oscillations observed in 
the torque signal. Finally, in section V we summarize the conclusions that our data make 
necessary. 

I. EXPERIMENTAL DETAILS 

Isothermal magnetoresistance sweeps have been performed using standard four probe 
techniques and a low frequency resistance bridge, on two single crystals at the National 
High Magnetic Field Laboratory, using the 33 T resistive magnet for the first and the 45 
T hybrid magnet for the second sample. The samples were pre-screened for In inclusions 
(by checking for superconducting signal in the resistivity at the T c of In) and mounted on 
a rotating sample holder inside the mixing chamber of a top loading dilution refrigerator. 
The sample 5*1 was measured with a current of 300 fiA up to 33 T at a fixed angle of 14° 
with respect to [001]. These data were later complemented by measurements up to 45 T 
on sample S2 with varying angle and a current of 1 mA. For both samples, the current is 
applied perpendicular to [001]. The field was rotated from [001] to a direction parallel to the 
ab-plane, but perpendicular to the applied current. The samples were allowed to thermalize 
for about 1 hr and the temperature of the mixing chamber was recorded at T (51) and 
12 T (S2) before each field sweep. Simultaneous torque magnetization measurements have 
been performed using an AC capacitance bridge on a sample cut from the first crystal (SI) 
and mounted with G.E. varnish on a Cu-Be cantilever. 

II. PHASE DIAGRAM 

The high field phase diagram has been previously established by magnetization 23 and 
specific heat^ measurements. Figures [T^i and [lb show the derivative of the torque signal as 
a function of magnetic field, H, for temperatures ranging between 0.04 and 0.92 K at fixed 
angle (—8°) and, in Fig. [lb, for angles (G) between H and [001] ranging from —1° to 29° at 
our base temperature (0.04 K). The torque signal is proportional to the vector product of 
magnetization with magnetic field and is shown for O = —8° and T = 0.04 K as an inset (see 
Fig. [lb). The torque signal is reported in Fig. [I] in arbitrary units since the cantilever used 
was uncalibrated. The feature in the torque signal, indicated by an arrow in figure [lb , or 



6 




20 25 30 35 



H(T) 

FIG. 1: Torque signal near the metamagnetic critical field, (a) Derivative of Torque signal vs. 
applied magnetic field, H, in Celrlns in the field range 20-35 T at an angle of —8° for temperatures 
ranging between 0.04 and 0.92 K. The peak corresponds to the metamagnetic transition, (b) Torque 
signal at 0.04 K (in arbitrary units) vs. H between 20 and 45 T for G = —8°. The metamagnetic 
field H m is indicated by an arrow, (c) Derivative of the Torque signal vs. H for field orientations 
ranging between —1° and 29° at the base temperature of 0.04 K. The data are shifted vertically 
for clarity. 

equivalently the peak in the derivative, Fig. [T^i and [Tfc, marks the metamagnetic transition 
field H m . H m is seen to shift to higher fields with increasing temperatures and angles. On 
the H — T phase diagram, shown in Fig. |2j the transition determined from torque data (at 
= —8° from [001]) extends to lower temperatures the previously reported transition line 
determined from the specific heat^. The transition field H m increases nearly quadratically 
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H(T) 

FIG. 2: Magnetic field phase diagram in Celrlns for i/||[001]. The metamagnetic transition H m is 
determined from torque magnetometry (•) and from specific heat data of Ref.— (*). The Fermi 
Liquid crossover temperature is determined from our previous specific heat data^ (o). The resis- 
tivity minimum for the two samples 51 (x) and S2 (+) are also shown. Inset: Angular dependence 
of the metamagnetic field H m determined from magnetization (•) and from resistivity (o) at 0.04 
K. The angle (0) is defined between the magnetic field and [001]. The dashed line correspond to 
a fit to H m = H m (0) rc^Q) . 

in temperature. This dependence of H m (T) is also observed in the Laves compounds and 
is consistent with Ginzburg-Landau theory^. The angular dependence of H m , shown in the 
inset of Fig. El can be fit to a co ^ e ^ dependence, suggesting that only the component of 
the field parallel to [001] induces the transition. In fact no transition is observed up to 60 
T for fields along the ab-plane^. The transition is most likely second order in the phase 
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space of temperature, field and angle that we explored, with no sharpening of magnetization 
anomaly into a first-order like discontinuity upon cooling, in contrast to the previous claim 
of Ref.- 25 . Although we cannot exclude a weakly first order transition, the torque data do 
not show any evidence for a critical end-point in Celrln 5 , as is the case for S^R^O^. 
The transition can be extrapolated to T = around H m (T = 0) = 28 T. Thus the phase 
diagram suggests the possibility of a metamagnetic QCP for H || [001]. Also shown in the 
phase diagram is the crossover line from Fermi Liquid to non-Fermi Liquid regime obtained 
in our earlier specific heat (C) investigation (up to 17 T) as the temperature below which 
j; becomes constant. 

III. MAGNETORESISTANCE 

The transition observed in M(H) is also apparent in the resistivity data shown in Fig. |3j 
The magnetoresistance (MR) in Celrln 5 is positive at low T, except around H m . Indeed, p 
goes through a broad maximum and decreases above H m for fields oriented close to [001] (0 < 
12°) (Fig. [3b). With a further increase of H, p increases steeply in the high field state. The 
field at which p is maximum closely follows the same temperature and angular dependence 
as H m determined from magnetization, and is also shown in the inset of Fig. |2j Note 
that a crossover from positive to negative MR has been reported in the related compound 
CeCoIn^, however, its relation to a possible metamagnetic transition is not well established. 
With increasing B the peak shifts to higher fields and becomes less pronounced: at G = 20° 
a plateau is observed. A small change in the slope can still be resolved at O = 27° but no 
hint of a transition is observed at = 34.5°. 

The temperature dependence of the resistivity is constructed from the isothermal sweeps, 
Fig. EJ and is shown for a range of fields in Fig. H] where p is plotted as a function of T 2 . 
The resistivity is consistent with a T 2 dependence between 0.04 K and 0.92 K at fields 
larger than 36 T where the sample has entered the field-polarized paramagnetic state. This 
is as expected for the case of a Fermi liquid ground state. Fits to a T 2 dependence to 
our data are shown as solid lines in Fig. for fields between 37 and 45 T. The slope in 
P — Po + AflT 2 defines the Fermi Liquid coefficient Apt,-, with Afl ~ 0.972 pQ cm/K 2 
(at H = 45 T) corresponding to a specific heat coefficient of 7 ~ 0.54 J/K 2 mol (using the 
Kadowaki- Woods ratio obtained for CeCoIn 5 in Ref.—), which is close to the experimental 



9 



7 

E 6 
o 

a 

CL 

5 
4 

20 25 30 35 40 45 

H(T) 



6 

E 
o 

CS 5 

=1 

Cl 

4 

3 

20 25 30 35 40 45 

H(T) 

FIG. 3: Magnetoresistance in Celrlns. Resistivity p vs applied magnetic field H in Celrlns for the 
sample 52 measured up to 45 T (a) for indicated temperatures at = —8° (b) for various angles 6 
between the field and the [001] orientation at the base temperature of 40 mK. Insets: (c) Angular 
dependence of p at T = 40 mK and H = 20 T in sample 52. (d) Resistivity vs applied magnetic 
field for the sample SI, measured up to 33 T, for indicated temperatures, at fixed angle = 14°. 

value found in Celrln 5 at H = 45 T (7 ~ 0.4 J/K 2 mol in Ref.— ), and slightly suppressed 
compared to 7 ~ 0.75 J/K 2 mol at H — 1 T 35 , as expected. This estimate suggests that 
the field-polarized state is still a heavy fermion metal, a conclusion also supported by the 
analysis of the cyclotron effective masses (see below). 

In contrast, for fields just above H m there is a strong deviation from a T 2 dependence of p 
at the lowest temperatures. Here the resistivity has a temperature dependence weaker than 
quadratic, as evidenced by the downward curvature in the figure for sample ^2 (Fig. Hb) be- 
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FIG. 4: Temperature dependence of resistivity in Celrlns for a field applied at = —8° from [001]. 
Resistivity, p, vs temperature, T, squared (a),(b) for the sample 52 between 12 and 45 T in the 
T-range 0.04 — 0.92 K and (c) for the sample 51 between 20 and 32.5 T in the T-range 0.03 — 0.87 
K. Solid lines in (a) are linear fits to p vs. T 2 . 

tween 33 and 35 T, and for sample SI between 30 and 32.5 T (Fig. Hb). This non quadratic 
T-dependence of p occurs very close to the critical field, H m , for both SI and S2. In a 
previous publication, where we measured a different Celrlns sample, we found that p was 
not well described by a Fermi liquid-like T 2 dependence for temperatures down to 50 mK 
and magnetic fields of 12, 15, and 17 T oriented parallel to the [001] direction whereas a 
T 2 dependence was found with the field oriented perpendicular to [001] 35 . This stresses the 
importance of magnetoresistance effects in the high field regime for the T— dependence of p, 
which can make the interpretation of non-Fermi Liquid behavior problematic. However, the 
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anomalous temperature dependence of p in Celrlns is clearly associated with the metamag- 
netic transition since the quadratic-in-T behavior is recovered for fields well above H m . Note 
that the occurrence of non Fermi-liquid like behavior in the resistivity is somewhat at odds 
with the simultaneous observation of quantum oscillations. This suggests that the inelastic 
(electron-electron) scattering responsible for the former becomes negligible compared to the 
elastic (disorder) scattering at very low T. This situation is not unique to Celrln 5 : CeCoIn 5 
also exhibits non-Fermi Liquid behavior concomitant to quantum oscillations close to its 
superconducting upper critical field2&~— Our observation of a non-quadratic temperature 
dependence of p near H m is in direct conflict with the observation of a resistivity well de- 
scribed by a T 2 dependence across the metamagnetic transition in CeRu 2 Si2^, suggesting a 
separate mechanism for IEM in these two heavy fermion systems. 

For fields below H m , a new contribution to p appears with a negative dp/dT below ~ 
0.5 K, as seen in Fig. HI This low-T upturn in p becomes more pronounced in the vicinity 
of the metamagnetic field, and is rapidly suppressed above it. The temperature T min where 
p(T) has a minimum is field dependent and defines the dome-shaped region in the H — T 
phase diagram, as shown in Fig. [2j We note that there is no anomaly associated with 
T min in the torque signal. This behavior has been observed in two different crystals at two 
different field orientations close to [001], and so it is a robust and reproducible feature of 
the metamagnetic transition in Celrln 5 . 



An upturn in p in the presence of high magnetic field is also observed at low— T in 
CeCoIng^ 1 ^ and UPt^. This contribution becomes more pronounced as the field is 
increased. It has been so far associated with the cyclotron magnetoresistance in the high 
field limit (u c t ~1), since it is only observed in the transverse geometry. An upturn is 
indeed possible in the high field limit of a compensated metal with p(H,T) ~ (oq+at 2 )* ' 
The upturn in Celrlns is saturating at low T and can be reasonably well described by this 
form. However, its suppression above the metamagnetic field is at odds with this scenario. 
This cannot be accounted for by the negative contribution to MR compensating the positive 
one, since MR is only negative close to H m and the upturn in p is not recovered when the 
MR becomes positive above 36 T. Moreover, the dome-shaped onset temperature T min , as 
shown in Fig. [2J clearly does not track the crossover to the high field regime (oj c t ~1), 
in striking contrast to CeCoIn^ 1 ^. To summarize, although MR effects in the high field 
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regime may account for the upturn in p(T), the suppression of this feature in the field 
polarized state, in the absence of Fermi Surface change(see below), suggests that a different 
mechanism is involved in Celrln 5 . 

In Sr 3 Ru207, a similar, but more pronounced, resistivity anomaly has been observed in the 
cleanest samples^ and has been tentatively assigned to the condensation of a nematic phase 
surrounding the QCP 21 . This is considered as a strong possibility in Sr 3 Ru207 since there 
are first order-like boundaries within which the resistivity anomaly occurs, making the idea 
of a phase transition appealing. The most likely candidate for such a domain-wall scattering 
scenario is Condon-like domains 2 ^, and the first order nature of the transitions makes the 
inhomogeneous state hypothesis a sensible one. Despite similar residual resistivities, the 
transition at H m in Celrlns is not as sharp as the one observed in S^R^Oy. Although we 
cannot rule out the possibility of a common mechanism, the absence of a sharp first order 
boundary make the domain wall scenario unlikely in Celrlns. 

An alternative explanation for the low-T resistivity anomaly near H m in Celrlns relies on 
the scattering of quasiparticles from quantum magnetic fluctuations leading to a crossover 
into a diffusive transport regime. This idea is consistent with the prediction that any devi- 
ation from Fermi Liquid behavior in presence of disorder will lead to an infinite resistivity 
at T = 41 . Moreover, logarithmic corrections to the conductivity, due to enhanced im- 
purity scattering in the presence of strong quantum fluctuations, have been predicted for 
a ferromagnetic 4 ^" 4 ^, and more recently for a metamagnetic 44 , quantum critical point. The 
idea is that near the QCP, the diverging magnetic coherence length (with decreasing T) 
eventually becomes larger than the mean free path (£) (or equivalently, the correlation time 
becomes larger than the time between successive collisions) and the system effectively be- 
haves as if it were impurity scattering limited. This is similar to what was previously reported 
in the itinerant ferromagnet Fei-^Co^Si 4 ^. The crossover from ballistic to disordered regime 
is expected to occur near T = ~ 0.62 K in Celrln 5 which is of the order of T min . Here 
we have used the Drude model along with the measured linear temperature coefficient of the 
specific heat (70 = 0.75 J mol -1 K -2 defined at H = 1 T), along with estimate of the electron 
density, 2.9 electron per Ce ion based on the Hall coefficient at 300 K— , n = 17.98 x 10 27 
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m 3 , in order to estimate the scattering time : 

r = ^- = ™ (1) 

ne 2 P (|)2/3(^)2 n 4/3 e 2 p ' ^ ) 

with the result that r ~ 1.2 10~ n s. However, as we lack convincing evidence of a QCP at 
H m in Celrln 5 , this hypothesis must be viewed as only a possibility. 

In summary, we believe several of the suggested mechanisms for the origin of the MR 
can be effectively ruled out, and one of them, quantum corrections in a diffusive transport 
regime, depends heavily on the existence of a QCP at H m , which we do not have conclusive 
evidence for at this time. In addition, the presence of non-quadractic T— dependence in 
p, which is clearly related with the metamagnetic transition in Celrln 5 , distinguishes it 
from the prototypical heavy fermion metamagnet CeRu2Si2- We note that in the latter, 
a metamagnetic QCP has been recently ruled out based on lower temperature, higher 
precision, resistivity measurements 3 -. The low— T upturn as well as the non-quadratic 
behavior in p raise important questions and present a challenge for understanding the 
unusual behavior that we observe in Celrlns. 



IV. DE-HAAS- VAN- ALPHEN EFFECT 

The study of the Fermi surface in heavy fermion metals via quantum oscillations of 
magnetization (dHvA effect) was pioneered by Taillefer et. al- 8 . These oscillations result 
from the discontinuity in the Fermi-Dirac distribution at T = and are well understood 
in simple metals in the framework of Lifshitz-Kosevich (LK) theory^. The observation of 
dHvA oscillations in heavy fermion metals is an important experimental milestone showing 
that these strongly correlated compounds can simply be described in terms of heavy quasi- 
particles, Landau quasiparticles with cyclotron masses much larger than band masses. The 
Onsager relation dictates that the dHvA frequencies are proportional to the cross-sectional 
area(s) of the Fermi surface perpendicular to the applied field. Therefore, by mapping the 
angular dependence of the frequency one can determine the Fermi Surface topography. 
This has been done in the 1-1-5 compounds, with the observed quasi-2D a (electron) and 
P (hole) FS sheets in close agreement with band structure calculations 3 ^^'. These results 
indicate that both CeCoIns and Celrlns have itinerant /-electrons at low temperatures, 
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FIG. 5: De Haas- van Alphen oscillations in Celrlns. Torque signal (after subtraction of a smooth 
background) vs inverse magnetic field H^ 1 at = —1° (a) and 29° (b). The corresponding field 
range is 25.6 to 40 T. The dashed line indicates the metamagnetic transition field H m =27.5 T and 
33.4 T for the two angles. 

whereas CeRhIn 5 has localized /-electrons. Moreover, a band-by-band investigation of the 
mass enhancement is possible with the dHvA technique, complementing the Fermi surface 
averaged effective mass given by the specific heat coefficient. By determining the angular 
dependence of the cyclotron mass it has been possible to identify "hot spots" in the parent 
cubic compound Celn^. 



Figure [5] shows the torque signal in Celrln 5 (in arbitrary units) as a function of inverse 
magnetic field H~ l , in the field range 25.6-40 T, at the base temperature 0.04 K, for two 
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FIG. 6: Angle dependence of the de Haas- van Alphen frequencies in Celrlns. (a) Frequency vs 
Angle in the low field (black symbols, H < H m ) and high field states (red symbols, H > H m ) 
corresponding to dHvA branches a\ (A), 02 (v)j a 3 (*) an d fa (•)• Zero angle corresponds to 
H || [001]. The fa orbit (at 10.4 kT) is only observed for = 14° and is not shown in this figure, 
(b) Fast Fourier Transform (FFT) spectra of the torque signal vs frequency in the low field (dashed 
lines) and high field (solid lines) states for field oriented at 29°, 22°, 14°, 7° and -1° from [001] (top 
to bottom). 

extreme angles —1° and 29°. The background field dependence (shown in Fig. db) was 
determined from a 5 th order polynomial fit to the torque signal and then subtracted from 
the data to better reveal the oscillations in the torque signal which are periodic in H^ 1 
corresponding to the dHvA effect. The dHvA frequencies are determined from the Fast 
Fourier Transform (FFT) spectra of the torque signal at each angle, shown in Fig. |6]d. The 
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FFT is performed on two field intervals of equal length in H~ l for H < H m and H > H m , 
thus keeping the same frequency resolution to allow comparison. Fig. [6^, compares these 
dHvA frequencies below and above the metamagnetic transition on an angular interval 
of ±30° around [001]. The observed frequencies correspond to extreme orbits on the a 
and (3 sheets of the Fermi surface, in agreement with previous reports (up to 16 T) in 
both their value and their angular dependence. The lower frequencies^ 1 ^, associated with 
the 3D electron pockets, are not easily detected in torque magnetometry since the torque 
signal is proportional to the angular derivative of the frequency (see Ref.—, p. 87). Overall, 
the same branches are observed both below and above H m , with only a slight shift in 
frequency across the metamagnetic transition. Moreover, within our frequency resolution 
(set by the field range of the FFT), we do not see any Zeeman splitting of the a nor 
peaks. This is generally the case in metals, as the Zeeman energy is small compared to the 
Fermi energyii^. These results rule out a major Fermi Surface reconstruction associated 
with the a and (3 Fermi surface sheets accompanying the metamagnetic transition in Celrlns. 

In the LK theory^, the effective cyclotron mass is determined from the temperature 
dependence of the amplitude of the dHvA oscillations. The dHvA amplitudes decrease with 
increasing temperature due to the thermal broadening of the Landau levels crossing the 
Fermi surface. We have determined the dHvA amplitudes by performing a FFT on intervals 
of equal length in inverse magnetic field if™ 1 centered around a field value ranging between 
20 and 33 T for the first data set (0 = 14°) and between 20 and 45 T for the second data 
set (0 = —8°). Each interval contains a fixed number of periods of a given dHvA frequency: 
29 periods of a 3 and (3 2 at —8° (45 T data set) and 23, 35, and 37 periods of a 3 , (3 2 and f3 1 
at 14° (33 T data set). The integer number of periods in a given field interval ensures no 
artificial broadening of the peaks with the periodic boundary conditions used in the FFT. 
The temperature dependence of the dHvA amplitudes is then fit to the LK formula in order 
to determine the effective masses for each branch at a given field (see Ref.—, p. 60): 

xT 

A = A 

sinn(xl ) 

a — 

with x = — where A is the dHvA amplitude, A Q the T = amplitude, m* the effective 
mass, m e the bare electron mass, and a = 2-** = 14.69 T/K. The data and fits shown 
in Figs. |7(a) and W(b) are from FFT spectra on a field range below H m (0) = 28 T. The 
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FIG. 7: Determination of effective masses from dHvA in Celrlns. (a) and (b) Normalized dHvA 
Amplitude, A, vs reduced temperature, T re( i = n ' , for orientations of = —8° and = 14° 
between H and [001], for the branches as(*), fai*), Pi(<>)- The amplitudes shown are from the FFT 
spectra on a small field range below H m . The solid lines are fits to data with the Lifshitz-Kosevich 
(LK) formula (see text) with Aq and m* being adjustable parameters. The amplitudes shown are 
normalized with respect to Aq, the T = value, (c) and (d) Effective mass vs magnetic field for 
G = —8° and = 14° determined from the LK fits. The dashed lines correspond to the IEM 
transition field H m . The vertical and horizontal error bars represent the uncertainty in the fitting 
parameter and the field range of the FFT respectively. The FFT is performed on intervals of equal 
length in H~ l centered around the indicated fields, corresponding to 29 periods for each branch 
for the 45 T data (a) and 23, 35 and 37 periods for 03,^2 and f3i for the 33T data (b). 
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field evolution of m* across the metamagnetic transition is determined in this fashion for 
each dHvA branch and is shown in Figs. [7(c) and[7J(d). The cyclotron mass decreases with 
increasing field for the (3% orbit, and is rather field independent for the other two orbits, 
and P>2- in contrast to other metamagnetic systems^ - — , no divergence of m* is observed 
in any branch for fields near H m , within the finite field resolution of our procedure set by 
the field interval of the FFTs. For the branches a.3 and 02 the effective masses show no 
significant angular variation between = —8° (Fig. [7(c)) and O = 14° (Fig. [7(d)). At —8°, 
the fix orbit is not resolved, possibly due to its large mass. The effective masses we obtained 
are also in close agreement with the previously reported masses for magnetic fields below 
16 T in the [001] direction^. 

The evolution of the Fermi surface has been studied via dHvA effect in a number of 
heavy fermion compounds across metamagnetic transitions^ 1 ^ 1 ^'^.. These studies have 
shown that the IEM is in general accompanied by drastic changes both in the effective 
mass and in the dHvA frequencies. This frequency change has been ascribed to /-electron 
localization in CeR^Si^ 1 ^ or to Zeeman splitting (Fermi surface polarization) in UPtjr^. 
In the light of these results, the absence of a substantial frequency change in Celrln 5 is 
striking, despite the fact that the field scale of H m is comparable to or larger than those of 
CeRu2Si2 and UPt3. This implies either a comparable or larger Zeeman energy in Celrln 5 . 
Recently, the IEM in Celrlns was attributed to a valence transition^ in which case we 
would also expect a change in the Fermi Surface volume. Our data appears to rule out this 
as a possibility. Thus, our data indicates that neither the Zeeman splitting, and associated 
proximity to a van Hove singularity, nor /-electron localization, or even valence fluctuations 
are the driving mechanisms for IEM in Celrlns. 

Despite the lack of drastic changes in neither the dHvA frequencies nor effective masses, 
we observe distinct changes in the dHvA signal near H m . Namely, the metamagnetic 
transition is accompanied by an overall damping of the dHvA oscillations. This can be 
directly observed on the background subtracted torque signal of Fig. [5j At = — 1° 
(Fig. [5^), only the lowest frequency 0:3 oscillations are observed above the transition field 
H m = 27. 5T, the 02 oscillations being completely suppressed in the field polarized state. At 
G = 29° (Fig. |5b), where 0:3 is the only frequency present in the entire field range, we see 
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FIG. 8: de Haas- van Alphen amplitude variation near the metamagnetic critical field in Celrlns. 
Contour plot of the FFT Amplitude vs. frequency and magnetic field, H, at T = 40 mK and (a) 
G = —1°, (b) 14°, and (c) 29°. The FFT spectra was obtained on a field interval centered around 
the indicated fields, all intervals being of equal length in H~ l , the length chosen to correspond to 
29 periods (panel a,c) and 23 periods (panel b) of 03. The a\ and 02 orbits are not well resolved 
due to the smaller field range of the FFT leading to broader peaks. The high frequency peak in 
(a) is the second harmonic of the 02 orbit. 

clearly that the amplitude of the oscillations does not increase monotonically with magnetic 
field but instead goes through a minimum around H m = 33. AT. The field evolution of the 
FFT spectra for three representative angles (—1°, 14° and 29°) is depicted on the contour 
plots of the dHvA amplitude vs magnetic field and frequency in Fig. [HJ The frequencies 
of the CV3, 0i and /3 2 orbits being all field independent, the corresponding contours stretch 
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vertically. As clearly seen in Fig. El the amplitude of the 02 orbit is completely suppressed 
at —1° in the field polarized state above H m . A similar trend is also seen for the 0\ peaks at 
14° (the data is only available up to 33 T at this angle). The 02 peak at 14° is not easily dis- 
tinguished from the neighboring cti, a<i peaks in the spectra due to the insufficient frequency 
resolution for small field ranges used in the FFT. In contrast, the amplitude of the a% peak 
is largest in the field polarized state after going through a minimum at a field just above H m . 

One can see this trend more clearly in the corresponding Dingle plots, shown for the same 
angles (—1°, 14° and 29°) as well as for several others in Fig. [9J Here the reduced dHvA 
amplitude as a function of inverse magnetic field H~ x is displayed. The reduced amplitude 
is defined as ln(A smfc ^P ), where A is the FFT amplitude of each orbit. According to the 
LK theory (see Ref.— , p. 66), the reduced amplitude should be linear in inverse field. As seen 
in Fig. [9^, this is not the case in Celrln 5 and the amplitudes of the 0i and 02 oscillations are 
strongly suppressed in the field polarized state for H > H m = 29 T at 14°. Moreover, the 
suppression of the 02 oscillations is not restricted to a particular angle but is systematically 
observed for the entire range where this orbit is resolved (—8° < < +14°), as shown in 
Fig. [9b. In contrast, the amplitude of the 0:3 oscillations go through a minimum above the 
transition but recover at higher fields, as shown in Fig. [9fc for various angles ranging from 
—28° to +29°. The maximum in the amplitudes of both 0:3 and 02 moves to higher fields 
(smaller H^ 1 ) with increasing angle, following the same trend as H m (see Figs. [9(b), [9(c)). 
Thus, the main result of our investigations is that there is a selective damping of the dHvA 
oscillations in the field-polarized state, due to the metamagnetic transition. 

In LK theory, the slope of the Dingle plot corresponds to the Dingle temperature defined 
as To = ^k^i w hh r _1 the scattering rate on each orbit describing the impurity broadening 
of the Landau levels at low-T— . For a 3 and 02 orbits, Dingle temperatures of 0.31 and 
0.27 K were obtained from linear fits to the reduced amplitudes in the low field state, for 
= —8°, using the effective masses determined independently from the LK fits (see above). 
These correspond to a mean free path of 112 nm for the a 3 orbit and 172 nm for the 02 
orbit which compares well to the value of 105 nm reported previously^. The scattering rate 
associated with such values of i is of the order of 3 THz, which is an order of magnitude 
less than the scattering time estimated from p . This is not surprising given the crude one 
band, one relaxation time approximation of the Drude formula as well as the differences 



21 



in the effective cyclotron masses and that estimated from the Sommerfeld model and 70 
(~ 260m e ) used in Eq. [TJ Nevertheless, the fact that the Dingle temperatures for both 
orbits are ~ 0.3 K which is of the order of both T and T min (see above) suggests that 
the transport is indeed in the diffusive regime with predominant disorder scattering. For 
0x, we could not determine Tjj since the Dingle plot is not linear in the inverse magnetic 
field, even in the low field state (see Fig. [9k). Thus the scattering rate for this orbit, as 
well as the effective mass (see above), appear to be somewhat field dependent, leading to a 
breakdown of the Dingle analysis. 

The observed decrease in the amplitude of the dHvA oscillations above H m for some of 
the observed branches is unexpected for at least two reasons. First, the spin fluctuations 
are suppressed in a field-polarized state which would indicate a smaller scattering rate, and 
thus larger oscillation amplitudes. Second, if one takes into account the magnetization in 
the dHvA effect, in fact the oscillations are in B' 1 rather than in H -1 , one would also 
expect a positive feedback with the larger internal field giving rise to a larger amplitude of 
the oscillations. A reduction of the amplitude for all of the orbits near H m may indicate an 
increased scattering rate due to electron scattering from magnetic fluctuations which are 
expected to be maximized at the phase transition. 

There are several reasons for a decreased dHvA amplitude at high magnetic fields that 
are understood from extensive dHvA experiments in simple metals. One such possibility 
for a decreased dHvA amplitude is the so-called magnetic interaction leading to phase 
smearing^, corresponding to the suppression of the overall signal due to electrons out-of- 
phase with one another in different parts of the sample. The underlying assumption is that 
the sample is either inhomogeneous in real space or that there is an inhomogeneous internal 
field distribution. In fact, if the magnetization is large then there will be an internal field 
inhomogeneity in the sample due to the shape dependent demagnetizing factor. Such 
an inhomogeneity would cause the dHvA amplitude to be smaller than when the field is 
homogeneous throughout the sample. This internal field correction is relevant when the 
magnetization of the sample is large. From the magnetization measured with a commercial 
SQUID at 1.8 K at an applied field of 7 T, we have extrapolated the volume magnetization 
in Celrln 5 to be about 326 G at 30 T (without taking into account the increase in M due 
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to the metamagnetic transition), which is 3 orders of magnitude smaller than the external 
field. Unless the sample itself is in an intrinsically inhomogeneous state, we are safe in 
assuming that B H is nearly homogeneous. To the best of our knowledge, no evidence 
for sample inhomogeneity is reported in the undoped, highly stoichiometric, high quality 
single crystals of Celrln 5 . Also it is not clear how this effect would lead to preferential 
damping of the /3 orbits, nor to the damping only in the field polarized state. Thus phase 
smearing is unlikely to be the cause of the observed damping. 

A second possibility for a non-linear Dingle plot is the magnetic breakdown for the ft 
orbits. Magnetic breakdown happens when the applied field is strong enough that the 
electrons tunnel from one piece of the Fermi Surface to another, leading to the observation 
of large dHvA frequencies^. Since we do not see any new frequencies appearing in the field 
polarized state, this is unlikely in our case. Finally, the dHvA amplitude can be suppressed 
due to the spin reduction factor, namely the destructive interference between the Zeeman 
split up- and down-spin electrons. Although we cannot resolve this splitting in the dHvA 
frequencies, such a polarization would lead to a decrease of the dHvA amplitudes. The 
spin reduction factor—, R s = cos( w v™* ), can be estimated for various orbits using the 
cyclotron effective masses determined above and with a Lande factor g = 2, yielding the 
same reduction factor for a and /3 orbits. Thus we can rule out Zeeman splitting as the 
cause of selective damping of dHvA oscillations above H m . 

Thus, one is lead to conclude that there has to be an additional damping mechanism 
beyond those typically considered (thermal and disorder broadening or Zeeman splitting) 
in the LK theory. Perhaps the most noteworthy aspect of our data is that the damping 
is selective: the heavier /3— orbits are more strongly damped than the lighter a— orbits, 
being totally suppressed in the field-polarized state. Although the cyclotron mass of all 
observed orbits does not change appreciably across the transition, a lighter heavy fermion 
ground state in the high field state would result due to this selective suppression of heavy 
quasiparticles from the heavier /3-sheet. It is also interesting to note that the /3 orbits 
are associated with the same hole band for which the smallest of the four hybridization 
gaps, namely A 1; has been assigned in optical conductivity measurements 33 . Although the 
interpretation of optical conductivity data may not be unique, it points to the same trend 
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in the evolution of this gap across the CeMIns family as the systematic changes found in 
other characteristic parameters such as 7 and T,M. 

The decreasing dHvA amplitudes for the f3— orbits above H m are also difficult to reconcile 
with the simultaneous drop in p. Assuming that p is determined by parallel conduction 
channels associated with each Landau orbit, and since the heavier /3— orbits (with smaller 
conductivity) are preferentially damped above H m , this would lead to an increased p, making 
necessary a simultaneous enhancement of the mean free path I of the light carriers in order 
to have a decreased p. Thus, the transport is most likely dominated by the lighter electrons 
on the a sheet that are strongly scattered near the transition but are otherwise intact in the 
field-polarized state. This is a possible explanation of the apparent paradox of a simultaneous 
decrease in both p and the dHvA amplitudes across the IEM transition. 

V. CONCLUSION 

In conclusion, our comparative study of p and dHvA oscillations is suggestive of an un- 
usual metamagnetic transition in Celrlns. The observed deviations from quadratic behavior 
in p vs T above as well as below H m , may well correspond to a non- Fermi Liquid regime 
extending over a wide range of magnetic field in the phase diagram. In the field-polarized 
state, above 36 T, a Fermi liquid regime is recovered. Moreover, an upturn is observed in 
the T— dependence of p that is somewhat different from what is reported in S^Ra^O^ 1 ^. 
The p(T) that we measure may correspond to quantum corrections that occur in a diffu- 
sively conducting strongly correlated metal. A detailed analysis of the dHvA oscillations 
shows that while the frequencies and effective masses are not significantly changed in the 
vicinity of the transition, the amplitude of some of the orbits are unexpectedly suppressed 
in the high field state. One possible way to reconcile the anomalous damping of the dHvA 
oscillations with the pronounced drop in the magnetoresistance is that the heavy holes are 
selectively damped across the IEM transition and that the transport is dominated by the 
lighter electrons. These results represent a challenge to our current understanding of IEM 
and call for further investigation as to what can cause such a change in the magnetization, 
the metamagnetic transition, beyond the mechanisms thus far explored, none of which can 
adequately account for the observed response in Celrlns. 
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FIG. 9: Dingle Temperature analysis in Celrlns. (a) Reduced amplitude ln(A stnh ^ x f^ ) vs. inverse 
field for @ = 14° for the orbits 03 (*), j3\ (A), 02 (•). For convenience, the corresponding 
field range is also indicated at the top of the figure, (b) Reduced amplitude vs inverse field for fa 
orbit at all angles where this frequency is resolved, angles ranging from —8° to 14° as indicated. 
The horizontal lines with arrow correspond to the noise level of the FFT spectra giving an upper 
limit for the amplitudes at the highest fields, (c) Reduced amplitude vs inverse field for 03 orbit 
at various angles. The data for all angles is up to 45 T, except for O = 14°, which is limited to 33 
T. The IEM transition field H m is indicated for each angle with a vertical arrow. 
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